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The vibrational modes of some single wall carbon nanotube (SWNT) intramolecular junctions 
(IMJs) have been calculated using the newest Brenner reactive empirical bond order (REBO) po- 
tential, based upon which their nonresonant Raman spectra have been further calculated using the 
empirical bond polarizability model. It is found that the Raman peaks induced by pentagon defects 
lie out of the G-band of the SWNTs, so the high-frequency part of the Raman spectra of the SWNT 
IMJs can be used to determine experimentally their detailed geometrical structures. Also, the in- 
tensity of the Raman spectra has a close relation with the number of pentagon defects in the SWNT 
IMJs. Following the Descartes- Euler Polyhedral Formula (DEPF), the number of heptagon defects 
in the SWNT IMJs can also be determined. The first-principle calculations are also performed, 
verifying the results obtained by the REBO potential. The G band width of the SWNT IMJ can 
reflect the length of its transition region between the pentagon and heptagon rings. 

PACS numbers: 63.22.+m, 78.30.Na, 61.46.Fg 



I. INTRODUCTION 

Single-wall carbon nanotube (SWNT) is a kind of 
nanoscale molecule obtained by wrapping a graphene 
sheet into a seamless cylinder. Because of its remark- 
able electronic properties^, SWNT is expected to play 
an important role in the nanoscale electronics, e.g., the 
field-effect transistor (FET) based logic devices 2 ^, and 
in applications of nanophotonics 4 . One of the most im- 
portant structures within these devices is the SWNT in- 
tramolecular junctions (IMJs), which are formed by in- 
troducing the pentagon-heptagon rings in them. Differ- 
ent chiral SWNTs can be connected through the IMJs 
and form the so-called metallic-metallic (M-M), metallic- 
semiconductor (M-S), or semiconductor-semiconductor 
(S-S) junctions. It is well-known that the electronic 
and optical properties of the SWNT IMJs have a 
close relationship with their geometrical or topologi- 
cal characteristics. A lot of theoretical studies on the 
electroniQ 5 ^ 6 ^ 7 ^ 8 ^ 9 ^ 10 ' 11 ! 12 ' 13 and optical^ properties of the 
SWNT IMJs have been carried out to get indirectly 
the detailed information on their topological structures. 
At the same time, many experiments have been per- 
formed to study the transport properties of the SWNT 
IMJs— i 16 i 17 i 18 . Recently, the direct Scanning Tunneling 
Microscope (STM) imaging and the measurements of the 
density of states (DOS) have provided the best direct 
evidences of the SWNT IMJsl2a2i. On the other hand, 
some efforts have been made to investigate the vibra- 
tional properties of the SWNT IMJs, e.g., the measure- 
ment of the radial breathing mode (RBM) by combin- 
ing inelastic electron tunneling spectroscopy and STM 
method 21 , and the Raman spectra by the confocal Ra- 
man spectral imaging 22 . Although the vibrational prop- 
erty obviously is very important for determination of the 
SWNT IMJs' structures, yet no theoretical study on their 
Raman spectra has been reported to date. 



Beyond the traditional force-constant method 23 , the 
state-of-the-art first-principles methods have already 
been used to calculate more accurately the phonon dis- 
persion curve a 24 i 25 i 26 i 27 i 28 , offering a unique technique to 
investigate the vibrational properties of the nanotubes 
or nanoropes. But for the very low symmetrical systems 
containing many atoms, such as SWNT IMJs, the first- 
principles method demands too intensive calculations. 
So, in this paper, we use the highly accurate empirical 
potential to calculate the vibrational properties of some 
SWNT IMJs. Here, the empirical potential is chosen to 
be the second-generation reactive empirical bond order 
(REBO) potential 29 , which is the newest version of the 
Tersoff-Brenner type potential, combining advantages of 
the two sets of parameters in the earlier version 30 . The 
first set of parameters (Brenner I) underestimates the 
isotropic elastic constants, while the second set of param- 
eters (Brenner II) overestimates the interatomic distance. 
The new REBO potential contains improved analytic 
functions and an extended database, offering a signifi- 
cantly better description of bond energies, bond lengths, 
and force constants for the hydrocarbon molecules, as 
well as elastic properties, interstitial defect energies, and 
surface energies for the diamond 2 ^. As a result, it should 
have a satisfying accuracy for describing the vibrational 
properties of the SWNT IMJs. 

After having obtained the vibrational modes of the 
SWNT IMJs, the nonresonant Raman intensity is fur- 
ther calculated by the empirical bond polarizability 
mode l 28 ! 31 . And for their Raman spectra, we will pay 
more attention to those characteristic Raman peaks 
caused by pentagon or heptagon defects because they are 
most important for investigating the local structures of 
the SWNT IMJs. 

In the following section, the REBO potential is first in- 
troduced and its ability to describe the vibrational prop- 
erties of SWNTs is discussed. Then, the vibrational 
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modes of some SWNT IMJs are calculated using the 
REBO potential. It is found that the characteristic Ra- 
man peaks caused by the pentagon defects lie out of the 
highest Raman peak of the SWNTs, which so can be re- 
garded as an indicator of the pentagon defects. In the 
later part of the paper, we will show that the heights of 
these peaks have some relationships with the number of 
the pentagon defects. On the contrary, the heptagon de- 
fects can only cause weak localized states below the G 
band of the SWNTs, which can be hardly observed in 
experiments. Furthermore, the junction length is shown 
to be able to influence the G band width of SWNT IMJ. 
In other words, the longer length of the transition re- 
gion between the pentagon and heptagon rings means 
the wider distribution of the G band. The first-principles 
calculations are carried out on the graphene sheet with 
Stone- Wales defect on it in order to further verify our re- 
sults quantitatively. Finally, this paper ends with some 
concluding remarks. 



II. NUMERICAL METHODS 



TABLE I: The Raman frequencies of (10, 10), (17, 0) and (5, 
5) SWNTs. The unit is cm -1 . Three different methods, i.e., 
the REBO potential, tight-binding method and LDA method, 
are used to obtain the dynamic matrix. The TB results are 
taken from Ref. l32l . The Raman mode around 900 cm -1 is 
E 2g in armchair tube and E\ g in zigzag tube, so in this table, 
it is labeled as E 2 (i) g - 
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Firstly, the phonon dispersion relationship and vibra- 
tional DOSs of armchair (10, 10) and zigzag (17, 0) tubes 
are calculated to demonstrate the accuracy of the REBO 
potential. Obtained results are shown in Figure 1. Cor- 
respondingly, the frequencies of the Raman active modes 
are presented in Table I. For comparison, the Raman fre- 
quencies calculated by both the tight-binding (TB) and 
the first-principles methods are also given. Here, the TB 
results are taken from Ref. 132 . 

The first-principles calculations are made as follows. A 
supercell geometry was adopted so that the SWNTs are 
aligned in a hexagonal array with nearby SWNT cen- 
ter distance of 25 A, which is found to be larger enough 
to prevent the tube-tube interactions. One unit cell of 
SWNT is sufficient to calculate accurately the phonon 




FIG. 1: (color online) The phonon dispersion relationship and 
vibrational density of states of armchair (10, 10) and zigzag 
(17, 0) SWNTs, calculated by REBO potential. Left part is 
for (10, 10) tube, and right one is for (17, 0) tube. 

modes at the T point. The k points sampling in the recip- 
rocal space is a uniform grid along the nanotube axis with 
the maximum spacing between k points being 0.03 A -1 
and the Gaussian smearing width is 0.03 eV. A plane- 
wave cutoff of 400 eV is used to obtain reliable results. 
After structure relaxation on both of the lattice constant 
along the tube axis and the atomic positions, the opti- 
mal structure is obtained when the residual forces acting 
on all the atoms were less than 0.01 eV/A. Our ah initio 
calculations were performed using highly accurate pro- 
jected augmented wave (PAW) method 33 , implemented 
in the Vienna ah initio simulation package (VASP) 34 , 
which is based on the density-functional theory in the 
local-density approximation (LDA). 

From Fig. 1 and Table I, one can conclude that the 
REBO potential can give reasonable lower Raman fre- 
quencies. As for the higher Raman modes, the REBO po- 
tential underestimates the slope of dispersion relationship 
(or the group velocity) and overestimates the frequen- 
cies. But its results agree well with those obtained by the 
TB method even quantitatively. In a word, the REBO 
method can be used to calculate the Raman modes of the 
SWNTs, especially for the tubes with larger radii. 

Next, we will introduce the construction method of the 
SWNT IMJs used in this work. 

The surface of SWNT and SWNT IMJs obeys the 
Descartes-Euler Polyhedral Formula (DEPF), i.e., for the 
open-ending SWNT or SWNT IMJ, the following equa- 
tion should be obeyed, 

V-E + F = 2{l-G). (1) 

Here, V, E, and F are the number of vertices, edges 
and faces, respectively, and x — 2 (1 — G) is the Euler 
characteristic. G is the genus, which is always equal to 
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FIG. 2: (color online) (10, 10)-(17, 0) SWNT IMJ. The o; 
timized structure is shown in the left panel and the period 
structure is shown in right panel. Only one pair of pentago: 
heptagon rings is used to connect the different tubes. 



1 for SWNT and SWNT IMJ with only one hole. F< 
perfect SWNTs, E = §V, F = \ , G = 1, and Eq. . 
is valid. If we further let N$ and Nj denote the number 
of pentagon and heptagon rings, respectively, then Eq. 1 
means that N 5 - ^7 = for any SWNT IMJ. In other 
words, the number of pentagon rings should equal to the 
number of heptagon rings. 

In Ref. 35 , a universal algorithm is carried out to con- 
nect arbitrary two SWNTs with one pair of pentagon- 
heptagon rings. Obviously, the structures generated by 
this method have the minimum defect energy. So, in this 
work, we also use this method to generate the SWNT 
IMJ structures. 



III. NUMERICAL RESULTS AND 
DISCUSSIONS 

Now, we firstly calculate the vibrational spectrum of a 
SWNT IMJ constructed by the (10, 10) and (17, 0) tubes 
with almost the same radii, whose structures are shown in 
Fig. 2. In order to eliminate the effect of dangling bonds, 
we further connect two same (10, 10)-(17, 0) SWNT IMJs 
to build a periodic structure shown in the right panel of 
Fig. 2, which is studied in this work. Then a structural 
optimization is performed to prevent the appearance of 
imaginary frequency, i.e., the soft mode. 

In Fig. 3, the nonresonant Raman spectra of some 
structures are presented. Because the low frequency R 
band and the high frequency G band are mostly inter- 
ested in experiment, only those two frequency ranges, 
i.e., 0-400 cm -1 and 1400-2100 cm -1 , are shown in Fig. 
3. 

Comparing Fig. 3a and 3c, it is clearly seen that many 
small peaks appear in the low frequency part. Some of 
them come from the original SWNT components, and 
their frequency numbers are affixed by letters of (i) or 
(c), indicating, respectively, the corresponding vibration 
mode comes from the in-phase or counter-phase combina- 
tion of the Raman modes in two SWNTs. For example, 
the peak at 25.2 cm -1 roots in the in-phase combina- 
tion of the E2 g modes of (10, 10) and (17, 0) tubes, but 
the peak at 34.3 cm -1 comes from their counter-phase 
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FIG. 3: (color online) The nonresonant Raman spectrum of 
the (10, 10)-(17, 0) SWNT IMJ. (a) Low frequency part, (b) 
High frequency part, (c) and (d) show the low and high fre- 
quency parts of the Raman spectra of (10, 10) and (17, 0) 
tubes, respectively. In Fig. 3(a), the letters of (i) and (c) 
indicate that the modes come from the in-phase and counter- 
phase combinations of the original Raman modes in the per- 
fect SWNTs, respectively. 




FIG. 4: (Color online) The atomic vibrations of the Raman 
active modes at frequencies of 25.2 and 34.3 cm -1 ), which 
come from the combinations of the basic Raman active modes 
of the (10, 10) and (17, 0) SWNTs. 



combination. Their atomic vibrations have been shown 
in Fig. 4, from which it is clearly seen that the atomic 
motions on the (10, 10) and (17, 0) tubes are in the same 
direction for 25.2 cm -1 mode, whereas they are in the 
reverse directions for 34.3 cm -1 mode. 

Some other peaks in the low frequency part, e.g., the 
peaks at 14.5, 20.6, 60.0, 65.5, 85.7 and 108.4 cm" 1 , can 
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FIG. 5: (Color online) The atomic vibrations of the Raman 
active modes at frequencies of 65.5, 85.7 and 108.4 cm -1 ), 
which come from the atomic vibrations of the junction part. 



not find their counterparts in the Raman spectra of the 
original SWNTs. Some of them come from the symmetry 
breaking of the perfect nanotube, and some others come 
from the vibration of the connecting part of the IMJ. For 
example, the atomic vibrations of 65.5, 85.7 and 108.4 
cm -1 modes are shown in Fig. 5. Obviously, the atomic 
motions near the defect parts are larger than those on 
the tube. But, most probably, it is practically difficult to 
observe these peaks because the most intensive low fre- 
quency peak found in experiments is the R band caused 
by RBM mode. 

As for the high frequency part, it can be found from 
Fig. 3b and 3d that many peaks appear in the G band 
(1600-1800 cm -1 ). A superposition of these peaks forms 
a widely distributed G band, which could not be regarded 
as a simple superposition of those G bands from the orig- 




FIG. 6: (Color online) The atomic vibrations of the Ra- 
man active modes beyond 1800 cm -1 for the (10, 10)- (17, 
0) SWNT IMJ. 



inal SWNTs. This phenomenon has been observed in the 
experimental result of Ref. [2l|, making also the detailed 
analysis on the G band of SWNT IMJ very complicated 
and hardly to be performed. 

On the other hand, in the range out of G band, i.e., 
in a range of about 1800-2050 cm -1 , some other Raman 
peaks can be found, whose modes are much simpler than 
the G band, and so may offer some useful information 
about the topological structures of SWNT IMJs. 

In order to understand the origin of these modes, their 
atomic vibrations are presented in Fig. 6. Obviously, 
these modes have a common character that their vibra- 
tional amplitudes are the largest on the pentagon ring, 
and become smaller and smaller on the other atoms of 
the tube surface when the atoms lie from the pentagon 
ring farther and farther. In other words, these modes 
are defect states caused by the pentagon ring. So, the 
pentagon can be regarded as a "topological impurity" 



in the perfect SWNT. Here, it is found that the bond 
length on the hexagon rings is about 1.42 A, while the 
mean bond lengths on the pentagon and heptagon rings 
are about 1.37 and 1.45 A, respectively. This fact means 
that the force constants on the pentagon (heptagon) de- 
fects should be larger (smaller) relative to those on the 
hexagon rings. Also, the distance between the pentagon 
and heptagon defects is quite large so that their inter- 
action can be neglected. So, we can simply use the one 
dimensional (ID) single impurity model in the classical 
solid state theory 36 to discuss qualitatively the vibra- 
tional property of here's (10, 10)-(17, 0) SWNT IMJ, 
and find that the frequency of the "topological impurity" 
state could be written as 
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Here, / and fo are the force constant of the "topolog- 
ical impurity" and the ID perfect lattice, respectively, 
with its maximum vibrational frequency of ljm(^m ~ 
4/o/M). And M is the atom mass. Eq. 2 is valid when 
/ > /oj which evidently corresponds to the situation of 
the pentagon defects. From Eq. 2, one can find that the 
defect frequency uo is higher than the maximum frequency 
of the perfect latticec^M, which is also well consistent with 
our numerical result. Furthermore, if we use the techni- 
cal term in the traditional mass-impurity model in the 
classical solid state theory^ 6 -, one may consider the pen- 
tagon defect to be the "light-mass impurity". But the 
frequency of the localized defect state can not be sim- 
ply obtained from Eq. 2 because it is difficult to make 
a direct link between the simple model and the realistic 
structure. 

Since the pentagon defect can cause localized defect 
state, then how about the heptagon defect? After search- 
ing for the whole Raman spectrum, we find that it is very 
hard to induce the localized states on the heptagon de- 
fect. The atomic vibrations of its two strongest defect 
states are shown in Fig. 7, which is found to be much 
weaker than those in pentagon defect states. Correspond- 
ingly, the Raman intensities of these modes are much 
weaker too. This phenomenon can be understood as fol- 
lows. Similar to the case of pentagon defect, the heptagon 
defect can be regarded as the "heavy impurity" in the tra- 
ditional mass-impurity model, which mainly causes the 
phase shift of the vibrational modes. And the frequencies 
of the localized modes caused by the "heavy impurity" 
are all lower than the maximum vibrational frequency 
of the perfect lattice, which is also consistent with our 
numerical result. 

To further validate our finding about the localized 
states, the first-principles calculation is made to study 
the vibrational properties of the pentagon-heptagon de- 
fects. But as indicated, the first-principles calculation 
on the SWNT IMJs is too expensive, and so is made on 
another substitute structure, i.e., a graphene sheet with 
a Stone- Wales (SW) defect on it. The SW defect is two 




FIG. 7: (Color online) The atomic vibrations of the localized 
states on the heptagon defect of the (10, 10)-(17, 0) SWNT 
IMJ. Here only the strongest two states are shown. 



fc> =1 505 cm" 1 




a =1726 cm" 1 



=1820 cm -1 

w 



FIG. 8: (Color online) Some localized states caused by pen- 
tagon defects on a graphene. 



pairs of pentagon-heptagon defects got together, and can 
be generated by rotating one bond in the graphene sheet 
by 90 ° . Now, the K-points sampling in the reciprocal 
space is a uniform grid in the graphene sheet and the 
maximum spacing between k points is still 0.03 A -1 . 

Firstly, the nonzero vibrational modes in a perfect 
graphene are calculated by the first-principles method, 
which is found to be 1591.38, 1591.14, and 896.24 cm" 1 , 
and similar to the results in Ref. l27|. Then, we have per- 
formed the calculation on a periodical (6x6) graphene 
unit cell with a SW defect in its center. The distance be- 
tween nearby graphene sheets is fixed to be 15 A, which 
is large enough to neglect the interaction between the 
sheets. Structural optimization is also made to prevent 
the soft modes. In such a structure, the distance between 
the SW defects is 14.5 A, which is also large enough to 
neglect the interaction between neighbor SW defects. 

The obtained results for three kinds of localized states 
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FIG. 9: (Color online) Some localized states caused by hep- 
tagon defects on a graphene. 
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FIG. 10: (Color online) Some states localized on the whole 
SW defects of a graphene. 



are shown in Fig. 8, 9 and 10. The localized states caused 
by pentagon defects are presented in Fig. 8, from which 
one can find their frequencies are very large, and those 
modes with larger frequencies tend to be more localized 
on the pentagon rings. For the localized states caused 
by heptagon defects shown in Fig. 9, the modes with 
larger frequencies also tend to be more localized on the 
heptagon rings, but their largest frequency is about 1356 
cm -1 , which is still lower than the G band frequency 
of graphene. All these phenomena agree qualitatively 
with our REBO results and the discussions based on the 
simple model. 

The localized modes in Fig. 10 can not be simply con- 
sidered to be localized on which defects. In fact, they 
distribute on the whole SW defect, which means they 
are caused by the whole SW defect. On the other hand, 
their appearance also indicates that the distance between 
pentagon and heptagon rings has a great inference on the 
behavior of localized states. 

Quantitatively, the highest frequency of the localized 
states on the pentagon rings in the SWNT IMJs should 
be higher than 1820 cm -1 due to the curvature effect of 
nanotubes. In fact, in Ref. [3?], the authors studied the 
vibrational property of the small diameter (3, 3) SWNT 




FIG. 11: (Color online) Two different structures of the (5, 5)- 
(10, 10) SWNT IMJs, viewed from different directions, (a), 
(b) and (c) are the high symmetry structures, (d), (e) and 
(f) are the low symmetry structures. 



with SW defect and found a local defect vibration mode 
with its frequency of 1962 cm -1 due to large curvature 
effect of the (3, 3) SWNT, which could be served as a 
fingerprint of the SW defect in carbon nanotubes. From 
our foregoing discussion, it can be deduced that this lo- 
cal defect vibration mode is just the highest frequency 
localized mode mainly caused by the pentagon defects. 

Now, we want to study the relationship of the local- 
ized mode frequency with the defect positions, and so 
make calculations on two structures, i.e., two (5, 5)- (10, 
10) SWNT IMJs connected by different methods, whose 
structures are shown in Fig. 11. 

The symmetry group of the first SWNT IMJ (Fig. Ha- 
lle) is C$ V1 which is called as the HS (high symmetry) 
SWNT IMJ. The symmetry group of the second SWNT 
IMJ (Fig. lld-llf) is C s , which is called as the LS (low 
symmetry) SWNT IMJ. Beside their different symme- 
tries, there are also another two important differences 
between them, which are the junction length defined as 
the distance between the pentagon-heptagon defects, and 
the number of pentagon- heptagon defects. There are 
five pairs of pentagon-heptagon defects in the HS SWNT 
IMJ, and only one pair of pent agon- heptagon defects in 
the LS SWNT IMJ. 

Now we will pay our attention only to the high fre- 
quency part of the Raman spectra of these two SWNT 
IMJs, which have been shown in Fig. 12. It can be found 
from Fig. 12 that the G band of the LS SWNT IMJ 
extends into a wider region than that of the HS SWNT 
IMJ. This should root in the long transition region in the 
LS SWNT IMJ, in which the IMJ radius and the bond 
lengths change almost continuously. So, most probably, 
the SWNT IMJ with a longer junction region will have a 
wider G band, which could be observed experimentally. 

In addition, one can find that the intensities of the 
highest frequency localized modes in the two SWNT 
IMJs are very different. In other words, relative to the G 
band intensity, the highest frequency localized mode in 
the HS SWNT IMJ has a larger intensity than that in the 
LS SWNT IMJ, which is probably due to the fact that 
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FIG. 12: (Color online) The high frequency part of the Raman 
spectra of the two SWNT IMJs. The upper and lower panels 
are for the high and low symmetry SWNT IMJ, respectively. 

there exist five times more pentagon-heptagon defects in 
the HS SWNT IMJ than in the LS SWNT IMJ. So the in- 
tensity of the highest frequency localized mode could be 
used to determine the number of pentagon rings exper- 
imentally. Once the number of pentagon rings is deter- 
mined, the number of heptagon rings can be determined 
too by Eq. 1. Generally speaking, nice reference sam- 
ples with a 'known' number of pentagons will be greatly 
helpful to determine the accurate number of defects in 
the IMJ, which is probably a technically difficult at the 
present time. However, with improving gradually the 
experimental tools, the nice sample could be manufac- 
tured in future. In a word, we think it is possible to use 
this method to determine the number of pentagon de- 
fects after we can set up a standard for it by measuring 
the nice reference sample with a " known" number of pen- 
tagons, e.g., a perfect SWNT with a cap at its end, which 
has six pentagons, or the SWNT with a Stone- Wales de- 
fect, which is composed of two pent agon- heptagon pairs 
next to each other. At least, probably this method could 
give an approximate number of pentagons in the SWNT 
IMJ. Finally, we also hope that it can be combined with 



other experimental methods, e.g., the optical responses 
to reach more practically this goal. 



IV. CONCLUSIONS 

In this work, the vibrational modes of some SWNT 
intramolecular junctions have been calculated by using 
the newest Tersoff-Brenner potential. Furthermore, us- 
ing the empirical bond polarizability model, the nonres- 
onant Raman spectra of these SWNT IMJs have been 
calculated. It is found that the highest frequency local- 
ized states caused by pentagon defects could be used to 
identify accurately detailed geometrical structure of the 
SWNT IMJ. That is because its frequency lies out of the 
G band of SWNTs and can be distinguished easily. And 
its intensity has a close relation with the number of pen- 
tagon defects. On the other hand, it is found that the 
heptagon defects can only cause weak localized states be- 
low the G band of SWNTs and can hardly be observed 
in experiments. All these results could be qualitatively 
explained by a simple "topological impurity" model. The 
first-principles calculations are also carried out to further 
prove our results quantitatively, and show that the char- 
acteristic frequency of the SW defects should lie between 
1820 and 1962 cm -1 , depending on the local curvature 
of the tube. Furthermore, the junction length is shown 
to have an important effect on the G band width of the 
SWNT IMJ, i.e., the shorter transition region will make 
the G band of SWNT IMJ to be narrower and simpler. 
Finally, the intensity of the Raman spectra has a close 
relation with the number of topological defects. 
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